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In a supernova core, magnetic fields cause a directional variation of the neutrino refractive index 
so that resonant flavor oscillations would lead to a deformation of the “neutrinosphere” for, say, r 
neutrinos. The associated anisotropic neutrino emission was proposed as a possible origin of the 
observed pulsar proper motions. We argue that this effect was vastly overestimated because the 
variation of the temperature over the deformed neutrinosphere is not an adequate measure for the 
anisotropy of neutrino emission. The neutrino flux is generated inside the neutron star core and 
is transported through the atmosphere at a constant luminosity, forcing the temperature gradient 
in the atmosphere to adjust to the inflow of energy from below. Therefore, no emission anisotropy 
is caused by a deformation of the neutrinosphere to lowest order. An estimate of the higher-order 
corrections must take into account the modified atmospheric temperature profile in response to the 
deformation of the neutrinosphere and the corresponding feedback on the core. We go through this 
exercise in the framework of a simplified model which can be solved analytically. 

PACS numbers: 97.60.Gb, 14.60.Pq, 98.70.Rz 


I. INTRODUCTION 

After the supernova collapse of a massive star, neu¬ 
trinos carry away about 99% of the gravitational bind¬ 
ing energy E\j of the nascent neutron star, taking with 
them a huge amount of momentum which is of the or¬ 
der 10"*^ g cm s“^ (i?b/3 X 10®^ erg). An anisotropy of 
the neutrino emission as small as 1% would suffice to ac¬ 
count for a neutron star recoil of about 300 km s~^ ||] 
and could thus explain the observed space velocities of 
most pulsars [^. However, even such a small asymmetry 
is difficult to explain. 

Pulsars tend to have strong magnetic fields, lead¬ 
ing to the speculation that B-fields could be a natural 
agent to cause asymmetric neutrino emission. For some 
time it appeared as if for realistic field strengths the in¬ 
duced polarization of the nucleon spins, together with 
the parity-violating properties of the neutrino-nucleon 
cross-sections, was enough to explain the observed pul¬ 
sar kicks Later it was recognized that this “cumu¬ 
lative parity violation effect” was in violation of funda¬ 
mental symmetries required of the Boltzmann collision 
equation; a correct derivation leads to a much reduced 
anisotropy [Q. 

This observation, together with the impressive recent 
evidence for neutrino oscillations, leads one to take se¬ 
riously another more indirect mechanism. The neutrino 
refractive index depends on the direction of the neutrino 
momentum relative to B. For suitable conditions, res¬ 
onant neutrino oscillations can occur between the neu¬ 
trinospheres of electron neutrinos and, say, r neutrinos, 
leading to a deformation of the effective Vr sphere [||, 
although the required conditions for large neutron star 


kicks may be rather extreme . The r neutrinos would 
thus be emitted from regions of varying effective tem¬ 
peratures, and thus, it was argued, would be emitted 
anisotropically. This idea was then taken up in sev¬ 
eral p^ers with modified neutrino oscillation scenar¬ 
ios 0-0. 

Unfortunately, however, this elegant scenario and its 
variations also appear to be fundamentally flawed in at 
least two serious ways. 

The first problem is caused by a common misunder¬ 
standing of the meaning of the “effective temperature” 
of the neutrino flux emerging from a supernova core. It 
is usually thought that the total energy carried away by 
neutrinos from a SN core is roughly equipartitioned be¬ 
tween the flavors, yet the heavy-flavor neutrinos (we usu¬ 
ally take Vt as an example) have stiffer spectra, i.e. their 
spectral temperatures tend to be much larger than those 
of Vg. Evidently, the neutrino luminosities are not given 
by the Stefan-Boltzmann law in a naive way—we will dis¬ 
cuss this issue in some detail in Sec. For the moment 
it suffices to observe that in a situation of exact flavor 
equipartition, a spectral swap of two flavors by oscilla¬ 
tions would not change the energy flux, except perhaps 
indirectly by a response of the thermal medium to the 
supposedly different spectra. 

The flavor equipartition of the energy flux need not 
be exact, and for the sake of argument we may contem¬ 
plate a situation where most of the energy is carried by Vg 
and i>g. If oscillations take place outside of the i^r sphere, 
the oscillated z/g’s could escape from different depths ac¬ 
cording to the B-field deformed resonance sphere, and 
thus with different temperatures. Even then one will not 
achieve a large flux asymmetry because it is not justified 
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to calculate the expected flux from the local gas temper¬ 
ature along the resonance surface. 

All of the neutrino spectra formation and oscillation 
physics of the present problem take place in the “atmo¬ 
sphere” of the protoneutron star, the outer region where 
the density drops quickly from core values around nu¬ 
clear density to effectively “zero”. The neutrino fluxes, 
however, are determined in the core of the neutron star. 
The atmosphere has virtually no heat capacity relative 
to the core. Therefore, after a short time, typically of 
the order of a few hundred milliseconds at most, which 
is very short compared with the Kelvin-Helmholtz neu¬ 
trino cooling time of the nascent neutron star, the neu¬ 
trino luminosity is governed by the core emission and the 
surface-near layers have reached a state where the tem¬ 
perature gradient ensures that all energy streaming up 
from below is carried outwards with a luminosity that is 
independent of the radial position. 

Therefore, the second serious problem of the oscillation 
kick scenario is that, to lowest order, a shift of the neutri- 
nosphere will leave the neutrino luminosity unchanged. 

A residual anisotropy effect obtains because the neu¬ 
trino flux is not strictly fixed by the core alone; it de¬ 
pends on the temperature at the core-atmosphere inter¬ 
face. This temperature, in turn, depends on the atmo¬ 
sphere so that there is an indirect influence of the atmo¬ 
spheric structure on the neutrino fluxes. More precisely, 
the neutrino flux determines the temperature gradient in 
the atmosphere, and the atmosphere influences the tem¬ 
perature at the core-atmosphere interface. Without a 
self-consistent treatment of this sort there is no pulsar 
kick at all, and the kick that one does obtain is a higher- 
order effect. 

In the following discussion we will elaborate our two 
arguments in more depth. In Sec. || we will explain the 
connection between the Stefan-Boltzmann law and the 
neutrino luminosities of a neutron star. We will stress 
the inadequacy of a simplistic application of the 
scaling of the luminosity when T is the spectral tempera¬ 
ture. In Sec. HI we will construct a simple self-consistent 
model in the so-called Eddington atmosphere approxima¬ 
tion Q . Our model leads to an estimate of the higher- 
order emission anisotropy from a changed neutrinosphere 
by direction-dependent resonant neutrino flavor conver¬ 
sions. Finally, Sec. 0 is given over to a discussion and 
summary of our hndings. 


II. NEUTRINO TRANSPORT IN NASCENT 
NEUTRON STARS 

A. Neutrino Eluxes and the Stefan-Boltzmann Law 

We begin our more detailed discussion with a descrip¬ 
tion of some crucial aspects of neutrino transport in 
nascent neutron stars. The picture thus developed will 
serve as background information for the analytical model 


of Sec. III. The most important insight to be presently 
explained is that the neutrino flux emerging from a super¬ 
nova core is not trivially given by the Stefan-Boltzmann 
law; the spectral temperature does not fix the flux, in 
contrast with true blackbody radiation [^. 

Lepton number is lost from the collapsed stellar core by 
the emission of electron neutrinos while energy is emit¬ 
ted in neutrinos and antineutrinos of all flavors. Elec¬ 
tron neutrinos are produced efficiently via the /3-process 
e~ + p ^ n + Vf. during the first second after collapse 
because of the high electron chemical potential, so that 
the deleptonization, in particular of the surface-near lay¬ 
ers, proceeds very fast. Most of the gravitational binding 
energy of the neutron star is radiated away after the col¬ 
lapsed stellar core has settled into the static, compact 
and hot protoneutron star when neutrinos and antineu¬ 
trinos of all flavors take up approximately the same share 
of the total energy and are emitted with very similar lu¬ 
minosities from the thermal bath of the core. 

The heat capacity and lepton number reservoir of the 
dense core are much larger than those of the less dense 
and much less massive atmosphere above. Roughly, core 
and atmosphere are discerned by the rather flat density 
gradient in the former, in contrast to the steep density 
decline in the latter. The density at the core-atmosphere 
interface is time-dependent and is typically between 10^^ 
and 10^^ g cm“^. Its small heat capacity and short neu¬ 
trino diffusion time imply that the atmosphere radiates 
away its binding energy in less than a few hundred mil¬ 
liseconds, to be compared with the neutrino diffusion 
timescale out of the core of a few seconds and the typical 
energy-loss timescale of several ten seconds. 

Thus, after a brief initial relaxation phase, the neutrino 
luminosity of the nascent neutron star is governed by the 
energy loss from the core; it reaches its surface value 
already below the core-atmosphere interface. Through¬ 
out the atmosphere, the luminosity is independent of the 
radial position if gravitational redshift is ignored. The 
temperature and density prohles in the atmospheric lay¬ 
ers adjust to the neutrino energy flux coming from inside 
to ensure its transport to the stellar surface under the 
constraint of hydrostatic equilibrium. This is equivalent 
to the situation in ordinary stars where the photon lu¬ 
minosity is produced in the nuclear burning zones while 
the stellar mantle and envelope adopt a structure in ac¬ 
cordance with the transport of this energy to the pho¬ 
tosphere. Evidently, the energy flux of neutrino or an¬ 
tineutrino species Vi cannot be given simply by the local 
gas temperature according to oc T'^. 

When thermal equilibrium between neutrinos and the 
stellar medium is assumed, the energy flux in the dif¬ 
fusion approximation can be expressed in terms of the 
atmospheric temperature gradient as 


Eh, = —D^ 


47r 

{hcf 
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Here, D^. is the diffusion coefficient, suitably averaged 
over the neutrino spectrum, while h, c and k are the 
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Planck constant, the speed of light, and Boltzmann’s con¬ 
stant, respectively. Further, 


^3 ) 
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1 -I- exp(a; - J ’ 


( 2 ) 


where rji,^ is the neutrino degeneracy parameter, i.e. the 
chemical potential divided by the temperature. In Eq. i) 
it was assumed that the neutrino degeneracy parameter 
is zero, =0, which then gives .^ 3 ( 0 ) ~ 3! = 6. This 
is always true for /x and r neutrinos and is also a good 
approximation for electron neutrinos and antineutrinos 
after the deleptonization of the atmosphere. 

Another expression for the neutrino energy flux can be 
obtained by relating it to the neutrino energy density, 

= ^(fcr)4-^3(0), (3) 


which yields 

F,^=c{fi)E,.,-^{kTf:Fs{0). (4) 


The factor (yx) E,^i denotes the average cosine of the angle 
of neutrino propagation relative to the radial direction. It 
is calculated from the neutrino phase-space distribution 
function (r, t, /i, e) according to 




f-i d/x yx /q°° de (r, t, yx, e) 

I-i fo°° de e^+Vi.i (r, t, yi, e) 


(5) 


where r is the radial position, t time, and e the neutrino 
energy. Choosing j = 0 gives us the “mean angle cosine” 
for the neutrino number flux, while j = I yields 

the corresponding quantity for the energy flux, (yx) E,vi ■ 
Comparing Eqs. ( 0 ) and gives 




dT 4 {X)e,., 
cT dr 3 hr 


( 6 ) 


where hr = {d\nT/dr)~^ is the temperature scale height 
and D^. = c(A) E,vil‘i was used for the diffusion constant, 
with the mean free path a suitable spectral average 
for the energy flux of neutrino species Vi. 

From Eq. (^) together with Eq. (^) one verifies that 
Fy. oc is not the whole story. Instead, the flux factor 
(yx) E,vi can be significantly different from the canonical 
value 1/4 which represents the Stefan-Boltzmann law. 
Put another way, the energy flux of “thermal” radiation 
is characterized by two parameters, the spectral temper¬ 
ature and the mean angle cosine which quantifies the de¬ 
viation from an isotropic phase-space occupation. 

Equation (||) reveals that (yx) E,vi depends on the posi¬ 
tion in the atmosphere because (A) E,vi becomes smaller 
for higher temperature and density. In the protoneutron 
star atmosphere the neutrino luminosity Ly = dxrr^Fh 
(for an individual neutrino type or for the sum of neu¬ 
trino and corresponding antineutrino) is fixed by the in¬ 
flow from the core region. The flux factor (yx) E,vi , on the 
other hand, increases with radius (decreasing tempera¬ 
ture) in accordance with Eq. (Q). 


B. Neutrino Spheres 

Another frequently misunderstood concept is that of a 
“neutrinosphere”. We stress that actually for each type 
of neutrino two different kinds of neutrinospheres are de¬ 
fined, the “energy sphere” with radius RE,vi and the 
“transport sphere” with radius Rt,yi- The latter is what 
many authors mean with the neutrinosphere, i.e. the sur¬ 
face of “last scattering” at optical depth 2/3 which emits 
the neutrino flux. The energy sphere is where energy¬ 
exchanging reactions freeze out while energy-conserving 
collisions may still be important. Of course, the concept 
of well-defined neutrinospheres or -surfaces is always a 
simplification of the real situation because the neutrino- 
matter interactions are strongly energy dependent. 

When RL,yi is the radius at which the luminosity Ly^ 
has reached its surface value, the three radii obey the 
relation RL,^i < RE,vi < Rt,vi, because the luminosity is 
fixed already deep inside the nascent neutron star. Be¬ 
tween RL,iii and RE,ui, the luminosity is constant while 
the spectral distribution of the flux still evolves due to 
neutrino absorption and reemission as well as energy¬ 
exchanging collisions. Between RL,ui and RE,vi the num¬ 
ber flux of a given neutrino flavor need not be conserved, 
whereas the lepton number flux (difference between neu¬ 
trinos and antineutrinos of a given flavor) is conserved, 
even for the electron flavor, because after a few hundred 
milliseconds the atmosphere is in a relaxed state and does 
not gain or lose lepton number on short timescales. 

Besides reactions in which energy is exchanged between 
neutrinos and the stellar medium, a significant fraction 
of the neutrino opacity of the stellar atmosphere is due to 
nearly iso-energetic neutrino-nucleon or neutrino-nucleus 
scatterings; for and Et this contribution in fact dom¬ 
inates. This has the consequence that even outside of 
Re,V i, where the flux spectrum is independent of radius, 
the neutrinos still undergo many scatterings and prop¬ 
agate outward by diffusion. Therefore, the local neu¬ 
trino distribution function is nearly isotropic, implying 
ih^j.vi 1- Hence, it is between RE,vi and Rt,vi, the 
region of a “scattering atmosphere”, where the naive 
Stefan-Boltzmann law is particularly poor at accounting 
for the neutrino luminosity [jl^ . 

The total opacity of electron (anti)neutrinos is domi¬ 
nated by /3-processes so that the distinction between the 
energy and transport sphere is often not crucial—for this 
flavor the concept of the neutrinosphere is crudely justi¬ 
fied. For the other flavors the distinction is crucial. 
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C. Pulsar Kicks by Oscillations? 

Our description of neutrino transport reveals that 
there is no simple relationship between the spectral tem¬ 
perature and luminosity of a given neutrino flavor [p^ . 
In numerical simulations with nonequilibrium transport 
description one finds approximately equal luminosities of 
neutrinos and antineutrinos of all flavors, but vastly dif¬ 
ferent spectral temperatures mi- However, the most 
elaborate numerical models have so far not fully taken 
into account several important energy-exchange channels 
between heavy-flavor neutrinos and the nuclear medium 
such as nucleon recoils, (inverse) nucleon-nucleon brems- 
strahlung, and collective as well as multiple-scattering 
effects 1^,1^; therefore, the spectral temperatures be¬ 
tween the flavors are likely far more similar than had 
been thought, and the difference between the energy and 
transport spheres may be less pronounced. 

In the limit of exact equipartition, neutrino oscillations 
along an aspherical resonance surface could not produce 
any pulsar recoil. One caveat is that the spectral swap 
between the flavors modifies their interaction rate with 
the medium so that there could be a small indirect effect. 
Even this possibility is diminished if the spectra are more 
similar than had been thought previously. 

Surely, it is not possible to calculate the pulsar recoil 
from a Stefan-Boltzmann flux, evaluated over the as¬ 
pherical resonance surface with its varying temperature. 
In typical simulations one finds an equipartition of the 
total energy to within a few percent. Any possible flux 
anisotropy caused by neutrino oscillations is therefore far 
smaller than had been assumed in Refs. [p[-p^. 


III. EDDINGTON ATMOSPHERE MODEL FOR 
FLUX ANISOTROPY 

A. Description of the Model 

In the preceding Section we have argued that in the 
limit of exact energy equipartition between the emitted 
neutrino flavors, an aspherical resonant oscillation sur¬ 
face could not cause a pulsar recoil, except perhaps by 
residual higher-order effects. But the equipartition need 
not be exact. This is especially true if the oscillations are 
into a sterile species which would not be emitted at all 
without oscillation effects Q. Even in this case it is dif¬ 
ficult to obtain a large recoil by oscillations because the 
magnitude of the anisotropy does not scale with the vari¬ 
ation of the gas temperature along the deformed emission 
sphere. 

An anisotropy of the neutrino emission can be estab¬ 
lished only in a much more indirect way. The aspheri¬ 
cal escape surface of the neutrino flux quickly leads to a 
perturbation of the initial configuration by producing an 
aspherical density and temperature profile of the neutron 
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FIG. 1. Schematic model of a nascent neutron star during 
the Kelvin-Helmholtz cooling phase. The core is treated as 
one zone, the atmosphere is described self-consistently by the 
Eddington atmosphere model. 

star atmosphere. This has a feed-back effect on the tem¬ 
perature at the core boundary and thus modifies the neu¬ 
trino flux. Therefore, one needs a self-consistent atmo¬ 
spheric model to estimate the neutrino flux anisotropy. 

To this end we subdivide the star into the “core” and 
the “atmosphere” as indicated in Fig. |^. For simplic¬ 
ity we take the “luminosity sphere” Rl as the interface; 
outside the neutrino luminosity is fixed. The core is char¬ 
acterized by a central temperature Tc, the temperature 
Ti at the core-atmosphere interface, and a linear tem¬ 
perature gradient in between. Clearly, the neutrino lu¬ 
minosity is determined by the difference between Tc and 
Tl- In order to estimate the flux anisotropy we need 
to calculate a self-consistent value for Tl by matching a 
self-consistent atmospheric model to the interface. 

Our main task, therefore, is to construct a model for 
the atmosphere which is characterized by three parame¬ 
ters: its mass, the temperature at the bottom, and the 
neutrino flux L which enters from below and depends on 
the central temperature Tc- The primary input quantity 
that varies as a function of direction is the atmospheric 
mass which is given by the neutrino resonance surface. 


B. Neutrino Eddington Atmosphere 

We construct a self-consistent atmospheric model by 
virtue of the Eddington approximation; for neutrinos 
this was done in Ref. , a fundamental paper that we 
will closely follow. The Eddington atmosphere employs 
the assumption of plane-parallel geometry, i.e. the atmo¬ 
sphere is taken to be geometrically thin relative to the 
core size. Moreover, one uses the diffusion approxima¬ 
tion, neutrinos and stellar medium are taken to be in 
thermal equilibrium, and the neutrino and antineutrino 
degeneracy parameters are taken to vanish everywhere, 
implying that there is no lepton number flux through the 
atmosphere. One uses neutrino and antineutrino opac¬ 
ities which are equal and vary with the square of the 
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neutrino energy, Ky. = Ap. = Ao(e^/eQ) with Aq be¬ 
ing a constant of dimension cm^g“^ and Eq = const. 
This implies that the neutrino and antineutrino phase- 
space distribution functions are identical, = /p^, and 
can be written in terms of the Fermi-Dirac distribution 
as 


= + (7) 

Here, /i is the neutrino angle cosine and m = dx p{x) 
the column mass density (g cm“^) of the atmosphere 
measured from the surface inward (Fig. 

These assumptions are reasonably well fulfilled for elec¬ 
tron neutrinos Ve and antineutrinos in the region of 
interest, i.e. between the and energy sphere Re,v^ 
and those of Vg and Vg, Re,v^ ~ Re^v^^ which are very 
close to the neutron-star surface. On the other hand, the 
thermal coupling of muon and tau neutrinos to the stellar 
background ceases in the relevant regions so that we take 
the atmospheric structure to be determined by the elec¬ 
tron (anti)neutrinos alone. Their combined luminosity 
Ai/e + = 47rr^(F^^ -|- Fp^) is assumed to be constant 

and given by the inflow from the core. In our plane- 
parallel model this implies that both the area AttP and 
the combined energy flux Fe = + Fp^ are constant. 

In this model one can derive an expression for the tem¬ 
perature as a function of the mass coordinate m which is 
given by Eq. (27) of Ref. ||ll| as 


Fe = 


2P 


CCn 


9 Ao{hc)^ rric 


[{kTgf - {kTLf 


( 11 ) 


When one defines the column mass density of the at¬ 
mosphere between radius Re and the surface at Rs as 
rriL ^ and plugs Eq. (H) into Eq. (^), one 

ends up with 


(kn)^ = a [{kTgf - {kn)^] 
+ v'/3[(A:Te)2-(fcFi)2] 


where 


rriL 


( 12 ) 


P = 


[kTgf 


20 


2l7r2 AqTOc 3 AR(Tc)mc 


(13) 


The second expression for /3 was derived by using Eq. (||) 
with T — Tg. The product Ar(Fc) mg is a measure of the 
optical depth of the core for Eg and Vg with the Rosseland 
mean free path being computed for a neutrino spectrum 
with temperature Tg. Equation (0) can be solved for r|, 

{kTL? = - P 

±PP[P + A{l + a){kTgf^ . (14) 


{kTf{m) 


/ 9{hc)^ ^ 
\ 2tt^ c e § 


FE\mF 


fZOihcf 
\ Ttt^ c 



1/2 


( 8 ) 


Since a <C 1 and f3 {kTgY (see below), this result can 
be approximated to first order in a by 

kTL = VakTg. (15) 


For m = 0 this equation yields the surface temperature, 
(fcTs)^ = Fe 30{hc)^PlPc), as a function of the energy 
flux Fe- Using Eq. (^ and the Rosseland mean opacity, 


Ar = 




(9) 


(Eq. 31 in Ref. |^]), one can write the energy flux Fe in 
terms of the derivative of the neutrino and antineutrino 
energy density e = as 


Fe 


c de 
3Ar dm 


( 10 ) 


(Eq. 28 in Ref. [|il). 

Equations (||) and ( 0 ) can be used to relate the en¬ 
ergy flux coming from the core to the temperature Te 
at the core-atmosphere interface. Taking Eq. (||) for the 
neutrino energy density with 2 ^ 3 ( 0 ) = Ttt'^/GO |1^] and 
TOc as the column mass density of the core, and evaluat¬ 
ing Eq. (^ in the one-zone approximation for an average 
core temperature, = ^{Tg + Tp) where Tg is the tem¬ 
perature at the center, one finds 


This and Eqs. & and (H) yield 

For conditions representative of the phase where the 
nascent neutron star loses most of its binding energy by 
neutrino emission, the central core temperature is around 
kTg ft! 30-50 MeV and the temperature at the base of 
the atmosphere (near or somewhat inside the muon and 
tau neutrino energy sphere) is /cTt ft 10-15 MeV. Equa¬ 
tion ( |r^ thus implies a ft 1/10 as a typical number. 
The optical depth of the core for Vg and Dg is of the or¬ 
der Ar(Fc)toc ~ 10®, so that P ^ 10~^{kTgp. Using 
these numbers in Eq. ( p^ one gets a luminosity Le = 
AttR^Fe ~ 2 X 10®i (Fs/10km)2(A:Tc/50MeV)'‘erg s-i 
for Vg plus Vg. Assuming all neutrino and antineutrino 
flavors contribute equally, this corresponds to a total neu¬ 
trino luminosity of ~ 6 x 10®® erg s“®, in good agree¬ 
ment with detailed numerical models. 
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C. Emission Anisotropy 


(23) 


In order to estimate the anisotropy of the neutrino 
emission we take the deformed effective neutrinosphere 
to be what corresponds to the atmospheric surface Rs in 
the previous section. The column mass density of the 
atmosphere, i.e. between the core boundary Rl and Rs, 
is taken as 


rriL = rriLo + Stul cos (j) 


m+ + m_ 
2 



COS(j) 


(17) 


where cos cj) = (q • B)/q is the cosine of the angle of 
the neutrino momentum q relative to the direction of 
the magnetic field and m± are the atmospheric column 
densities that correspond to co^= ±1. With Eq. ( 0 ) 
inserted into Eq. (|^) and Eq. (|iq ) one obtains 

Fe(cos (j>) (X P {kTcY f 1 — cos (()^ . (18) 

\ rric nic ) 


The asymmetry in the third component of the neutrino 
momentum can now be estimated as 


Ag 1 d cos (j) cos p Fe (cos p) 
^ ® cos (j) Fe (cos p) 

1 SrriL 

18 TOc 


(19) 


This result is accurate to first order in ao = tolo/wc 1 
and we have assumed, as in Refs. |||] and Q, that only one 
neutrino species is responsible for the anisotropy which 
carries off about 1 /6 of the total energy. 

The mass difference <5 tol is connected with the radial 
deformation Sr of the surface of resonance, which is de¬ 
fined by r{(j)) = ro + Sr cos p, through 

SruL = ^{m+- m_) = poSr . (20) 


Here, po is the mean density at the surface of resonance. 
For the width Sr in dependence of the strength of the 
magnetic field B one finds |5|j|] 





-1 


3 eR 

2 ^^ 


( 21 ) 


where = \d\nne/dr\~^ is the scale height for 

changes of the electron number density near rg and 
pe ~ /ic( 3 ne/ 87 r)^/^ is the chemical potential of the elec¬ 
trons. With the definition 



and the density scale height hp = \d\r\p/dr\^^ near rg, 
Eqs. (p0| ) and ( |^ ) yield 


SrriL = Pohp 7 . 


Finally, with Eq. (^^ one ends up with 
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1 Pq hp 

T?- 

18 rric 


Taking pohp/mc < rriLc/mc ~ 1/10 and hn^/h^ 
leads to the numerical estimate 


< 


Aq 

— < - 0.0012 


20MeVY / B 

Pe ) Vl0l®G 


(24) 
1 § 

(25) 


This result is at least 10 times smaller than the aniso¬ 
tropy derived in Ref. [^. 

Therefore, the kick mechanism based on a deforma¬ 
tion of the effective neutrinosphere requires more than 
an order of magnitude larger magnetic fields than esti¬ 
mated in Refs. ii-i whose analysis already reduced the 
effect originally discussed by Kusenko and Segre [|| . For 
a neutrino emission anisotropy of 1 %, corresponding to 
a recoil velocity of the nascent neutron star of approxi¬ 
mately 300kms~^, one needs magnetic fields in excess of 
about 10^^ G near the stellar surface. 

The value of 7 in Eq. (1^) is sensitive to the electron 
chemical potential. In Ref. 7 was evaluated by using 
Ye = rie/nb « 0.1 for the electron fraction {ub is the 
baryon number density) at a density p ~ 10^^ gcm“^. A 
discussion of the uncertainties of this choice can be found 
in Appendix A where the structure of the protoneutron 
star atmosphere is self-consistently determined from a 
simple, analytical model. Typically, 7 decreases during 
the neutrino cooling of the nascent neutron star because 
the atmosphere becomes denser and more compact as the 
star cools and deleptonizes. This disfavors large Sr at 
intermediate and late times during the Kelvin-Helmholtz 
cooling when most of the gravitational binding energy 
of the neutron star is emitted in neutrinos. This makes 
large emission anisotropies even more unlikely. 


IV. SUMMARY AND DISCUSSION 

We have argued that the neutrino-oscillation scenar¬ 
ios for neutron star kicks i -0 suffer from two serious 
flaws; both problems are related to an incorrect picture of 
neutrino transport in the atmosphere of a protoneutron 
star. 

First, when the neutrino luminosity is equipartitioned 
between the flavors, no significant recoil can be pro¬ 
duced because only an indirect, higher-order effect re¬ 
mains which is associated with the spectral swap of two 
neutrino flavors. Spectral differences imply a change of 
the neutrino interaction with the stellar background and 
thus affect the neutrino transport from the core to the 
surface through a modified atmospheric temperature pro¬ 
file. While we cannot estimate the magnitude of the re¬ 
sulting small kick velocity, we are convinced that it is a 
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very small effect. Moreover, current supernova models 
overestimate the spectral differences between Vf. and Vt 
because neutrino interactions have not been taken into 
account which enhance the thermal coupling between Vj- 
and the stellar medium 

Second, when the luminosities are taken to be vastly 
different (as had effectively been assumed in previous pa¬ 
pers), again no effect obtains in zeroth order because a 
shift of the location of the neutrinosphere does not affect 
the neutrino luminosity. The latter is governed by the 
core emission, not by local processes in the atmosphere. 

However, the atmosphere adjusts to the oscillation- 
induced modihcation of its transport capabilities, causing 
a small, higher-order effect due to an altered temperature 
gradient in the core; in our treatment it was expressed as 
a change of the temperature at the core-atmosphere in¬ 
terface. In this sense the neutrino flux from the core hxes 
the atmospheric temperature profile and, conversely, the 
atmosphere determines the core emission. In contrast, in 
Refs. it had been assumed that the unperturbed 

atmospheric temperature is a measure of the neutrino 
luminosity by virtue of the Stefan-Boltzmann law. 

Both of our arguments imply a huge suppression of the 
pulsar recoils calculated in Refs. 1-0 , but a realistic 
quantitative estimate of the residual effects is not possi¬ 
ble with our simple analytic tools. A detailed numerical 
treatment would be extremely difficult, and the motiva¬ 
tion for such an effort is minimal because most likely 
one would conhrm what now looks like a non-effect. In 
any case, it is clear that the oscillation scenarios require 
much larger magnetic fields than had been contemplated 
in Refs. |5|-|l0f| and thus probably take one beyond what 
is astrophysically motivated. 

We find it disappointing that both the cumulative par¬ 
ity violation and the neutrino oscillation scenario, which 
seemed to work with reasonable magnetic field strengths, 
do not survive a self-consistent discussion. Of course, it 
remains possible that huge magnetic helds 10^® G) 
with an asymmetric distribution in the core of the pro¬ 
toneutron star cause sufficiently asymmetric neutrino 
opacities for a large neutrino rocket effect [|^. It is also 
possible that asymmetric neutrino emission has nothing 
to do with the pulsar kicks. Either way, it does not look 
as if the pulsar velocities can be attributed to neutrino 
oscillations within presently discussed scenarios. 

Note Added in Proof. —In a recent preprint [21], 
Kusenko and Segre criticize our analysis and afhrm their 
previous results [5]. Their main objection against our 
work is that allegedly we ignored neutrino absorption via 
charged-current interactions and assumed equal opacities 
for all neutrino flavors. However, opacity differences be¬ 
tween electron neutrinos and muon/tau neutrinos were, 
of course, included in our analysis. In our analytical 
model, the opacities determine the column density of the 
atmosphere between the core boundary and the neutron- 
star “surface”. The latter was taken to be the effective 
sphere of neutrino-matter decoupling which is located at 
different radii for the different neutrino flavors. 


In their new analysis [21], Kusenko & Segre estimate 
the neutron star kick associated with anisotropic reso¬ 
nant flavor conversions by considering the asymmetric 
absorption of electron neutrinos in the magnetized neu¬ 
tron star atmosphere. This approach is based on the 
same assumptions as their previous one [5] and there¬ 
fore it is not astonishing that their original estimate of 
the magnitude of the pulsar kick is confirmed. However, 
resonant flavor conversions in the neutron star atmo¬ 
sphere cannot cause a persistent emission or absorption 
anisotropy because the neutrino luminosity is determined 
by the flux from the core, and the atmosphere adjusts to 
the inflow from below within a time which is very short 
compared to the neutrino-cooling time of the nascent 
neutron star. Therefore, the absorption anisotropy cal¬ 
culated by Kusenko and Segre is a transient phenomenon 
until the enhanced absorption is balanced by the reemis¬ 
sion of neutrinos and the atmosphere has approached a 
new stationary state. Of course, if the magnetic held is 
initially present rather than being “switched on,” even 
this transient phenomenon will not occur—the unper¬ 
turbed atmosphere simply never exists. 

As stressed in the main text of our paper, the large 
kick velocities found by Kusenko and Segre are an arti¬ 
fact of using an unperturbed atmospheric model instead 
of a self-consistent one. Accepting that a young neutron 
star is well described by our core-atmosphere picture (and 
Kusenko and Segre do not seem to question this crucial 
premise of our work), our conclusion that there is no 
zeroth-order kick caused by neutrino oscillations in the 
atmosphere is rigorous and as such not subject to de¬ 
bate. As described in our paper, there will be a higher- 
order effect due to a modihcation of the temperature at 
the core-atmosphere interface caused by the asymmetry 
of the self-consistent atmospheric structure. The result¬ 
ing kick is much smaller than estimated by Kusenko and 
Segre. Our discussion, however, does not exclude that 
a larger, globally asymmetric neutrino emission may de¬ 
velop if the anisotropies are produced by effects in the 
dense inner core of the neutron star. 
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APPENDIX A: SIMPLE MODEL FOR 
PROTONEUTRON STAR ATMOSPHERE 

A crucial parameter for estimating 5r in Eq. j^ ) 
is the electron number density in the protoneutron 
star atmosphere between the energy spheres of and 
heavy-havor neutrinos and v-r- Making use of the 
fact that the degeneracy parameter of electron neutri¬ 
nos approaches zero, rji,^ = rje + Vp ~ Vn ^ Q, as the 
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protoneutron star atmosphere deleptonizes and the elec¬ 
tron chemical potential ^’e = kTrje decreases one 

can easily estimate Ue and the electron number fraction 
Ye = Ue/nb = (n-e- — Ue^^jub where nb = plrriu is the 
number density of baryons and rriu the atomic mass unit. 
With Ye = Yp, Yn = 1 - Yp, and rin - rjp ^ hi{Yn/Yp) 
for Boltzmann gases of neutrons n and protons p with 
mn ~ TUp, one finds as a very good approximation 


_ SirmuikT)^ Ve , 2 , 2 ^ 

{hcfp 3 ^ 


9.1 X 10 


-3 


kT Y 
5 MeV / g cin“^ 


-1 




inl 1 


n 2 


(Al) 


This equation shows that Ye decreases with increasing 
density p for a given temperature. 

The temperature as a function of column mass den¬ 
sity m is given by Eq. (S as = Am + (kTsY. In 

order to determine p{rn), we employ hydrostatic equi¬ 
librium which yields for the pressure P = gm + Ps as 
a function of m. Here Pg is the pressure at the sur¬ 
face and g = GM^s/Rs is the gravitational acceleration 
near the surface of the protoneutron star which is nearly 
constant in the thin, plane-parallel atmosphere. Since 
the layers between the energy spheres of electron neutri¬ 
nos and muon and tau neutrinos are dense (p ^ 10^^- 
10^^ g cm“^) and rather cool {kT 3-10 MeV), the 
pressure is dominated by baryons (see Fig. 8 of Ref. |^|) 
so that we can take P « kTp/mu to obtain 


p{m) 


{gm + Ps)mi, 
kT 


(A2) 


With Eqs. (B|) and (|A2|) one derives 


{kTf _ [Am+{kTgf\ 


mugm + pskTg 


(A3) 


From this relation and Eq. one can see that Ye first 
decreases only slightly, then increases with rising m, i.e. 
on the way inward into the atmosphere. Therefore, the 
minimum value of the electron fraction can be found very 
close to the electron neutrino energy sphere RE,ue ~ Rs- 
This confirms hn^/hp < 1 because p as well as Ye have 
a negative gradient interior to Re,^^- Typical conditions 
near Re,v^ at early times during the protoneutron star 
evolution are |ll5| , ^ kTg « 3 MeV and ps ~ 10^^ g cm“^ 
which yields Ye{Rs) ~ 0.078. At later times the tem¬ 
perature in the atmosphere drops due to cooling |]l^ , |l^ , 
and according to Eq. (A2) the density in the atmosphere 
must become higher and the density gradient steeper. 
Therefore, the electron neutrinosphere, which is assumed 
to be located at a certain value of the optical depth, 
moves to higher densities. Typical conditions then are 
kTg ~ 5 MeV and ps ~ 10^^ g cm“^ for which one gets 


Ye{Rs) ~ 0.050. Even later, one has kTg ~ 5 MeV and 
Ps « 5 X 10^^ g cm“^ which gives Ye{Rs) ~ 0.004. 

Although Ye{Rs) decreases as the protoneutron star 
cooling goes on, the number density ne{Rs) nevertheless 
increases because of the rising density ps. In the listed 
examples, ne{Rs) changes from 4.7 x 10^^ cm“^ through 
3.0 X 10^^ cm“^ to 1.2 x 10^® cm“^. From Eq. ( |^ ) we 
conclude that this disfavors large Sr at intermediate and 
late times during the Kelvin-Helmholtz neutrino cooling 
of the nascent neutron star. 


[1] S.E. Woosley, in: D.J. Helfand and J.-H. Huang (eds.), 
The Origin and Evolution of Neutron Stars (Kluwer, Dor¬ 
drecht, 1987), p. 255. S.E. Woosley and T.A. Weaver, in: 
D. Pines, R. Tamagaki and S. Tsuruta (eds.). The Struc¬ 
ture and Evolution of Neutron Stars (Addison-Wesley, 
Redwood City, CA, 1992), p. 235. 

[2] A.G. Lyne and D.R. Lorimer, Nature 369, 127 (1994). 
D.R. Lorimer, M. Bailes and P.A. Harrison, Mon. Not. 
R. Astron. Soc. 289, 592, (1997). J.M. Cordes and 
D.F. Chern off, Astrophys. J., in press (1998), astro- 


ph/9707308. 


[3] 


[4] 


[5] 

[ 6 ] 

[7] 

[ 8 ] 

[9] 

[ 10 ] 

[ 11 ] 


C.J. Horowitz and Gang Li, Ph ys. Rev. Lett. 80, 3694 
(1998). D. Lai and Y.-Z. Qian, astro-ph/9802344 (ver¬ 
sion 1), withdrawn; Astrophys. J. 495, L103 (1998); (E) 
ibid. 501, L155 (1998). 

A. Kusenko, G. Segre ar id A. Vilenkin, |astro-ph/9806205 . 


P. Arras and D. Lai, astro-ph/9806285, submitted to 


4872 


Phys. Rev. Lett. (1998). 

A. Kusenko and G. Segre, Phys. Rev. Lett. 77, 
(1996); ibid. 79, 2751 (1997). 

Y.-Z. Qian, Phys. Rev. Lett. 79, 2750 (1997). 

A. Kusenko and G. Segre, Phys. Lett. B 396, 197 (1997). 
E.Kh. Akhmedov, A. Lanza and D.W. Sciama, Phys. 
Rev. D 56, 6117 (1997). 

D. Grasso, H. Nunokawa and J.W.F. Valle, Phys. Rev. 
Lett. 81, 2 412 (1998). 

R. Horvat, hep-ph/980638C. 


P.J. Schinder and S.L. Shapiro, Astrophys. J. 259, 311 
(1982). 

[12] H.-Th. Janka, Astrop. Phys. 3, 377 (1995). 

[13] R.W. Mayle, J.R. Wilson and D.N. Schramm, Astrophys. 
J. 318, 288 (1987). R.W. Mayle, in: Supernovae, ed. 
by A.G. Petschek (Springer, New York, 1990). J.R. Wil¬ 
son and R.W. Mayle, 1989, in: The Nuclear Equation of 
State, Part A, ed. by W. Greiner and H. Stocker (Plenum 
Press, New York, 1989). S.W. Bruenn, Phys. Rev. Lett. 
59, 938 (1987). S.W. Bruenn, 1993, in: Nuclear Physics 
in the Universe, ed. by M.W. Guidry and M.R. Strayer 
(lOP, Bristol, 1993). H.-Th. Janka, in: Prontier Objects 
in Astrophysics and Particle Physics, ed. by F. Gio- 
vannelli and G. Mannocchi, (Societa Italiana di Fisica, 
Bologna, 1993). H.-Th. Janka, Ph.D. Thesis (Technische 
Universitat Miinchen, 1991), MPA-Report 587 (Max- 
Planck-Institut fur Astrophysik, Garching, 1991). 


8 























[14] H. Suzuki, Num. Astrophys. Japan 2, 267 (1991). H. 
Suzuki, in: Frontiers of Neutrino Astrophysies, Proc. of 
the Internal. Symp. on Neutrino Astrophys., Oct. 19-22, 
1992, Takayama/Kamioka, Japan, ed. by Y. Suzuki and 
K. Nakamura (Universal Academy Press, Tokyo, 1993). 

[15] H. Suzuki, in: Physics and Astrophysics of Neutrinos, ed. 
by M. Fukugita and A. Suzuki (Springer-Verlag, Tokyo, 
1994). 

[16] H.-T. Janka, W. Keil, G. Raffelt and D. Seckel, Phys. 
Rev. Lett. 76, 2621 (199 6). S. Hannestad an d G. Raffelt, 
Astrophys. J., in press, astro-ph/9711132. A. Burrows 


and R. Sawyer, Phys. Rev. G 58, 554 (1998). 

[17] S.A. Bludman and K.A. Van Riper, Astrophys. J. 224, 
631 (1978). 


[18] G.S. Bisnovatyi-Kogan, Astron. and Astrophys. Transac¬ 
tions 3, 287 (1993). E. Roulet, J. High Energy Phys. 1, 
13 (1998). D. Lai and Y.-Z. Qian, Astrophys. J. 505, 844 
(1998). 

[19] W. Keil and H.-Th. Janka, Astron. A strophys. 296, 145 
(1995). H.-Th. Janka and W. Keil, astro-ph/9709012, 
MPA-Report 1043 (1997). J.A. Pons, S. R eddy, M. 


Prakash, J.M. Lattimer and J.A. Miralles, astro-ph/ 


980704C 


[20] S.E.Woosley, J.R. Wilson and R. Mayle, Astrophys. J. 
302, 19 (1986). 


[21] A. Kusenko and G. Segre, astro-ph/9811144 


9 







